Dynamical electron transport through a nanoelectromechanical wire in a magnetic 
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We investigate dynamical transport properties of interacting electrons moving in a vibrating na- 
noelectromechanical wire in a magnetic field. We have built an exactly solvable model in which 
electric current and mechanical oscillation are treated fully quantum mechanically on an equal foot- 
ing. Quantum mechanically fluctuating Aharonov-Bohm phases obtained by the electrons cause 
nontrivial contribution to mechanical vibration and electrical conduction of the wire. We demon- 
strate our theory by calculating the admittance of the wire which are influenced by the multiple 
interplay between the mechanical and the electrical energy scales, magnetic field strength, and the 
electron-electron interaction. 
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Recent progresses in experimental techniques of nano 
structures have made it possible for researchers to in- 
vestigate electromechanical systems with extremely small 
sizes, [lj The size limit is now being pushed down to the 
vicinity of mechanically quantum regime and there are 
much effort to understand the effect of quantum fluctua- 
tions on electrical and mechanical properties — and their 
coupling — in the nanoelectromechanical systems. @, d, 
EL S S 0, S] Recently Shekhter and his coworkers in- 
vestigated an Aharonov-Bohm effect in quamtum me- 
chanically vibrating nanoelectromechanical wires. In 
the limit where non-interacting electrons tunnel between 
weakly connected leads through virtual processes, they 
used a perturbative analysis to find that the quantum 
interference has a strong effect on the electromechanical 
response of the wire. 

In this Letter, we propose an exactly solvable theo- 
retical model of strongly interacting electrons in a nano 
wire which quantum-mechanically oscillates in a mag- 
netic field. We study the dynamical electron transport 
beyond perturbative regime and investigate the funda- 
mental limit of both quantum mechanics and strong cor- 
relations. Whereas Ref. focuses on energy scales much 
less than the level spacing of the vibrational normal 
modes, our model is valid for a much broader range of 
energy scales, including energies much greater than the 
vibrational energy level spacing. Therefore, it captures 
various important and interesting physical consequencies, 
such as admittance resonances occuring near mechanical 
normal mode frequencies. We also consider the electron- 
electron interaction in the theory, which is known to play 
an important role in low dimensional systems such as 
carbon nanotubes. [13 Our theory is based on a Lut- 
tinger liquid model which was previously developed by 
the present authors. Q One of the important character- 
istics of it is that it treats both electrical and mechani- 
cal degrees of freedom quantum mechanically, completely 
on an equal footing. On the atomic level, the interplay 



between electrical charge and mechanical oscillation of- 
ten causes rich and interesting new physics to emerge 
- the most famous examples being the BCS theory of 
superconductivity[ll| and the SSH theory of conducting 
polymers. [l2j| This Letter deals with a mesoscopic ver- 
sion of such a case, in which one can even control the 
electromechanical coupling by tuning external parame- 
ters. 

FigureQ]schematically shows the experimental setup of 
our model. A one-dimensional wire — a single- wall car- 
bon nanotube, for example — of length L is suspended 
across a valley between two metallic gates. There are 
two important degrees of freedom for this model: the 
mechanical oscillation of the wire and the electric cur- 
rent through it. Mechanically, the wire may form stand- 
ing waves of which sound velocity is determined by its 
one-dimensional mass density and tension. When the 
amplitude of the oscillations is small, the mechanical mo- 
tion of the wire may be characterized by its displacement 
u(x) from the equilibrium position, where x is the one- 
dimensional position along the wire. By connecting the 
gates to an external ac voltage source, one may now drive 
an ac current through the wire. If there is a magnetic 
field B perpendicular to the wire, the wire is expected 
to feel the Lorentz force and move. This may in turn 
induce electromotive force in the wire as it cuts through 
the magnetic flux. This is the well-known effect of back- 
reaction. [3} We will discuss this effect in the regime where 
both the electron current and the mechanical oscillations 
must be treated quantum mechanically. 

The electronic excitations of our one-dimensional elec- 
tromechanical system is best described by the Tomonaga- 
Luttinger liquid theory. [l3l 14, [l5| In this theory, the 
electron-electron interaction cannot be treated perturba- 
tively, but the Hamiltonian may be rendered quadratic 
by a bosonization technique. This is the key to the exact 
solvability of our model. For simplicity, we will consider 
a spinless case here, but it may be easily generalized to 
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include spin. The electronic part of the Euclidean action 
is given by 1|| 17 1 
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where vf is the Fermi velocity, r = it is the imag- 
inary time, and [3 = 1/fcsT is the inverse tempera- 
ture. Here, 8(x) is a bosonic field related to the one- 
dimensional electronic density fluctuation Sn via Sn — 
(d8/dx)/2ir. We have replaced a system of finite-sized 
one-dimensional Tomonaga-Luttinger liquid wire con- 
nected to two-dimensional leads by an effective one- 
dimensional liquid with position dependent interaction 
parameter and velocity: [la] 



K(x) = 



K, if < x < L, 

1, if x < or x > L. 



(2) 



For a short range interaction of the form V(x — x') — 
VS(x—x'), we have K = 1/ yl + V/nvp. Inside the wire, 
the velocity of the acoustic plasmon is also renormalized 
to v = vf/K. In a more realistic interface, K(x) would 
change more smoothly, but it will not make qualitative 
changes to the results below. 

The mechanical part of the action is 
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where u{x) is the transverse displacement of the wire 
from the equilibrium position, p the one-dimensional 
mass density of the wire, and v s the sound velocity of 
the mechanical transverse waves. 

Finally, the two fields 9 and u are coupled via a mag- 
netic field B = Bz. We will choose a guage in such a 
way that the vector potential is given by A = —Byx. 
At a given position along the wire, one may simply re- 
place y with the transverse position of the wire u(x). The 
coupling part of the action is then given by 

Se-u = - I dr I dxJA x 
c Jo Jo 
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where — e is the electron charge and J = (—e/2ir)(d8/dt) 
is the elctric current. Since the total Action 
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is completely quadratic, this model is exactly solvable. 

Before going into any detailed calculations, we would 
like to first discuss the qualitative characteristics of this 
model. In the absence of magnetic field, Sg. u = 
and the electrical and mechanical degrees of freedom 



are completely decoupled. Then there are two inde- 
pendent branches of excitations. The normal modes 
of the mechanical degree of freedom u(x) are stand- 
ing waves u m (x) — u m0 sin(rmrx / L) , the energies of 
which are quantized as mA u with an integer m. Here, 
A u = irfiv s /L is the energy level spacing of the mechan- 
ical normal modes. On the other hand, the electrical de- 
gree of freedom 9{x) is infinitely extended and its energy 
is not quantized. However, due to the facts that the volt- 
age drop occurs within the length of the wire and K(x) 
changes at the boundaries, the energy scale Ag = nhv/ L 
may still manifest itself in some physical quantities as will 
be shown below. This in fact is the energy level spacing 
of the Tomonaga-Luttinger acoustic plasmon confined in 
a box of length L. 

If the magnetic field is now turned on, the two degrees 
of freedom become coupled to each other through the 
following two processes: (i) a current induces a Lorentz 
force on the wire and (ii) an oscillation causes electromo- 
tive force through magnetic induction. Therefore, stim- 
ulating one of the two degrees of freedom will induce 
excitations in another. Eventually, the effect will nega- 
tively feed back into the originally stimulated degree of 
freedom; this is the origin of the back-reaction. What 
characterizes the coupling strength is the magnetic back- 
reaction energy scale lob = By/vFe 2 /irhpc 2 . This is the 
size of the back-reaction gap in one of the two uncou- 
pled excitation branches. [3| For a small magnetic field, 
the coupling is weak and the electronic excitations are 
only slightly perturbed. The effect of the coupling will 
be most strong at each quantized energy level of the me- 
chanical oscillations, which, in the weak field limit, occurs 
at every A u . (In fact, even multiples of A u will have no 
effect due to the reason explained below.) As the mag- 
netic field increases, the coupling grows stronger and the 
two branches of excitations intricately merge together to 
form new sets of excitations. The dispersion relation for 
the new sets of excitations have been calculated for an 
infinitely long wire in Ref.II In the strong coupling limit 
(ojb 3> uj Ae/A u ), the dispersion relations are approxi- 
mately given by E(q) w v s vq 2 /lob and ^/v 2 q 2 + uj%. As 
we will see below, these will characterize the energy pro- 
file of finite wires, too. 

One of the important physical quantities that charac- 
terize the system is the current response to an external 
ac bias voltage. Let us consider the average current in- 
side the wire I = J Q J(x)dx. The external bias volt- 
age and the current will be assumed to take the form 
V cxt = V Q {oj)e ZUJt and I = I a (uj)e %UJt , respectively. Then, 
we may calculate the ratio of the complex amplitudes of 
the current and the external bias voltage in the linear 
response regime, i.e., Y(u) = lim Va ^)^ 1 (w)/V (u>), 
which is the inverse of the impedance and is simply called 
the admittance. For small biases, the linear response the- 
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ory dictates 
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Y(w) = -—. dx dx' dte luJt U(x,x',t) 7 (6) 
iTiujL- 1 J J J 

where 

U(x,x',t) = (J(x,t)J(x',0)) = (^) 2 (6{ Xl t)e{x' 

(7) 

is the current-current correlation function. This may be 
computed using the Euclidean action in Eq. ([5]) and the 
usual analytic continuation technique. 

The real part of the admittance is proportional to the 
power absorption spectrum. Our calculations show that 



Re Y(lo) 
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r\ = A u /A g = u a /«. 
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These quantities need to be evaluated numerically. 

Figure [2] shows the main result. Assuming that the 
wire is a carbon nanotube, we have taken the value 
K = 0.22 from Ref. [l(]. Since the mechanical sound 
velocity is usually small compared to the electron Fermi 
velocity, we have also assumed r\ = A u / Ag = v s /v = 0.1. 
For ujb — 0, the mechanical degree of freedom is decou- 
pled from the electronic one and does not contribute to 
the current measurement at all. Therefore, Ag is the only 
relevant energy scale and it determines the frequency 
scale over which KeY(u>) decays at small u>. It also deter- 
mines the period of weak modulation in the admittance, 
which is usually too small to see in the figures. Note 
that the admittance approaches the usual Landauer dc 
conductance e 2 /h, as to — > 0. 18[ This dc conductance is 
unaffected by magnetic field and stays unchanged for all 



values of ujb, simply because there is no back-reaction 
for a time- independent dc current. 

For a weak magnetic field (lob <C <-u/ri), regularly placed 
sharp resonance peaks start to appear. Their positions 
are determined by the condition T(uj) — 0. Note that 
they occur whenever the frequency uj is close to an odd 
integer multiple of A u /h. From the positions of the res- 
onances, we can deduce that the changes in the mechan- 
ical excitation energy levels are perturbatively small in 
the weak coupling limit. The widths of the peaks are 
usually very small, with Q factors sometimes reaching as 
large as ~ 10 7 . They increase with the magnetic field as 
B 2 , which may be easily understood as broadening ef- 
fect due to the electromechanical coupling. There is no 
resonance at even integer multiple of A u /Ti due to the 
following reason. The magnetic flux swept by the mth 
normal mode is B u m {x)dx 1 but this always vanishes 
for an even integer m because the areas above and below 
the wire cancel each other. 

As the magnetic field grows, the coupling becomes 
stronger. The positions and shapes of the resonance 
peaks change much and some peaks even get wiped 
out. In the limit of strong magnetic field's ^ ^l^i)-, 
there are well distinguishable sharp peaks at u « 
m 2 AgA u /hcuB- This is in good agreement with the 
dispersion relation for the lower-energy gapless branch 
of the strongly coupled system with q = mir/L.^ If 
frujB ^> Ag,A u , there are also broad peaks that develop 
near u ^ lob, which is a result of crossover between the 
weak and strong coupling regimes. 

Now let us discuss about the effect of the electron- 
electron interaction. From Eq. ((8]), it is easy to see 
that ReY(to) depends on K only through an overall 
coefficient of T(co). In general, Re Y(ui) decreases as 
the interaction becomes stronger, except right on res- 
onance where ReY"(ti;) becomes independent of the in- 
teraction strength because T(lo) = 0. Well away from 
resonances [r(w) 3> uj], the value of KeY(oo) is approxi- 
mately proportional to K 2 . 

All quantities discussed so far have no temperature de- 
pendence, because the action of our theory is completely 
quadratic. It will not be true if, for example, the con- 
tacts between the wire and the leads are not perfectly 
transmitting. In that case, the action will no longer be 
quadratic and there will be nonvanishing thermal effects 
in general. For partially transmitting contacts, we also 
have to consider the effect of charge quantization and 
Coulomb blockade, although carefully adjusting a back- 
gate voltage may lift Coulomb blockade. [l9j Although 
Ref. y studies an opposite limit of weakly coupled leads, 
the results agree with ours in the sense that the electron 
transmission is depressed by magnetic fields, because it 
is a direct consequence of the back-reaction. 

In summary, we have shown that the quantum mechan- 
ically coupled electronic and mechanical degrees of free- 
dom of a one-dimensional wire oscillating in a mangetic 
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field may be studied in an exactly solvable model. Calcu- 
lating the admittance, we have found that the interplay 
of three important energy scales, the oscillation energy 
level spacing A„, the Tomonaga-Luttingcr liquid energy 
level spacing Ag, and the magnetic energy lob may re- 
sult in rich and interesting consequencies such as sharp 
resonance peaks and crossover in terms of relative magni- 
tude of the frequency to the magnetic induction coupling 
strength. From the positions and shapes of the reso- 
nance peaks, we may extract such pieces of information 
as the magnetic field, electron-electron interaction, and 
the standing wave energy levels, which may be used in 
applications such as nano sensors. 
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FIG. 1: Schematic figure of the setup. A nanowire (thick 
solid line) is suspended between two metallic gates and os- 
cillates about its equilibrium position. It is influenced by a 
perpendicular magnetic field B. 




uj/A u w /A. 

FIG. 2: Left panel shows the density plot of the real part 
of the admittance ReF as a function of u)/A u and lob/A u 
where ojb = B^/vFe 2 /nhpc 2 , and A u = nhv s /L, and u) 
is the ac bias frequency. The parameters used here are 
r) — A u /Ae = Vs/v — 0.1 and K — 0.22. Lighter regions 
correspond to higher values. The curves in the right panel 
show cross sections of the density plot at several different 
magnetic fields; from bottom to top, ujb/A u = 0, 0.2, 9.8, 
and 26.5. Curves are offset for better viewing. 



